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Abstract 

The effect of topology on the thermodynamics of a gas of adjoint representation 
charges interacting via 1+1 dimensional SU(N) gauge fields is investigated. We demon- 
strate explicitly the existence of multiple vacua parameterized by the discrete superse- 
lection variable k = 1, • • • , N. In the low pressure limit, the k dependence of the adjoint 
gas equation of state is calculated and shown to be non-trivial. Conversely, in the limit 
of high system pressure, screening by the adjoint charges results in an equation of state 
independent of k. Additionally, the relation of this model to adjoint QCD at finite 
temperature in two dimensions and the limit ./V — > oo are discussed. 
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1 Introduction 



It has been known for a long time now that a gauge theory coupled to matter can have non- 
trivial vacuum structure. For our purposes this statement is taken to mean that in the theory 
there exists a family of isolated (vacuum) sectors of which the physical one must be chosen 
by a superselection rule. The classification of such '^-vacua' can be carried out in all cases by 
examining the topological structure of the effective gauge group which acts non-trivially on 
the matter content of the theory. Unfortunately, this classification has little to say about the 
consequences of the choice of vacuum on the physics of the system. Our purpose here is, in 
the case of two dimensional SU (N) gauge fields coupled to a gas of classical adjoint charges, 
to explicitly determine the thermodynamics of the system as function of vacuum sector. 

The SU (N) non-Abelian Coulomb gas with adjoint charges in two dimensions is a useful 
model for investigating the topological and symmetric properties of gauge theories in higher 
dimensions. In the limit of infinite N this model has a first order phase transition which 
can be interpreted as a deconfining one separating a phase with tightly bound fundamental 
degrees of freedom and a plasma-like phase. This behaviour is analogous to the situation in 
higher dimensional Yang-Mills theory and adjoint QCD where the transition is characterized 
by a breaking of the center symmetry of the gauge group and the Polyakov loop serves an 
effective order parameter. Unfortunately, for the case of finite iV which we will be studying 
here there is no such phase transition in two dimensions. However at finite N the vacuum 
structure due to the topological structure of the theory is clearly apparent. Moreover the 
model is explicitly solvable in the limits of high and low particle density and in these limits 
we will investigate the thermodynamics of the adjoint gas in each sector of the theory. 

The standard method of classifying the multiplicity of vacua in a particular gauge theory 
with adjoint matter depends on identifying the effective gauge group. Here, since gauge 
transformations operate by adjoint action on all fields, a transformation which lies in the 
center Z is trivial. Consequently the true gauge group is the quotient of the gauge group and 
its center. This quotient is multiply connected, as can be seen from the relation, 

TT^G/Z) = 7T (Z) = Z (1.1) 

where G is any semi-simple gauge group. Thus, Z gives a classification of gauge fields which 
are constrained to be flat connections at infinity. In that case 

lim A^x) = ig\x)V^g{x) (1.2) 

|x|— >oo 

where g(x) is a mapping of the circle at infinity to the gauge group GjZ. Since G/Z is a 
symmetry of the Hamiltonian, we expect that all physical states carry a representation of 
7Ti(G/Z). In the case where the center of the group is Abelian all of its irreducible repre- 
sentations are one dimensional and further, when Z ~ 7L ni x • • • x 7L n . , we are lead to a 
classification of all physical states in terms of j generators of Z, {zi, • ■ • , Zj}. If Z is a unitary 
realization of Z and \if) > is a physical state we have 

Z\i)>= e i{zi+ - +z ^\ij > (1.3) 
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The equivalence class of all states which have the same transformation properties under Z 
form an isolated (vacuum) sector of the theory which is typically called a 9 sector. This 
somewhat abstract motivation for the existence of multiple sectors in a theory is in complete 
agreement with more concrete constructions of 9 sectors in 1 + 1 dimensional Yang-Mills theory 
without matter content |], H, [|. The objective of the current paper is to develop a framework 
for extending this analysis to systems with arbitrary matter content with particular attention 
given to the case of matter in the adjoint representation of the gauge group. 

The thermodynamics of classical charges interacting via Abelian and non-Abelian electric 
forces in one spacial dimension has been considered previously [|, [|. Additionally the effect of 
multiple vacua in the 377(2) adjoint gas has been considered previously ||. We will consider 
constant pressure ensembles as these authors but we will use a different formalism to construct 
the partition function of the system that can be easily extended to configurations other than 
that of the open line. 

We begin in the next section with a short description of our methods for constructing the 
model of 1 + 1 dimensional non-Abelian Coulomb gas. Restricting ourselves to the case of the 
adjoint charges, we proceed with an analysis of the low density /pressure limit of the model. 
Here using group theoretic techniques the explicit dependence of the equation of state of the 
adjoint gas on the vacuum parameter k is established. Converting to the Fourier domain we 
find the high density /pressure limit of the model is equivalent to solving a system of coupled 
quantum oscillators. In this limit the equation of state is shown to be independent of k. We 
conclude with a discussion of the limit of large rank groups, iV — > oo, and the connection of 
our model to adjoint two dimensional QCD. 

2 The Classical Non-Abelian Coulomb Gas 

The physical system which we will be investigating is that of a collection of static charges 
restricted to lie on a line, interacting with each other via 1+1 dimensional non-Abelian electric 
forces. Since we are interested in the thermodynamics of this system we will consider time 
to be a compactified coordinate with period equal to the inverse temperature (3 = l/T. 
Consequently the topology of interest is that of a cylinder. In this section we demonstrate 
the construction of a 1 + 1-D Coulomb gas of non-Abelian charges using group theoretic ideas 
introduced by Migdal and developed later by Rusakov ||. From this point of view the 
propagator for colour electric charge along a spacial distance L with boundary holonomies 
given by the unitary group elements Ui and U2 is, 

K[Ui, U 2 ; A] = J2xr(Ui) e~ A X r(U 2 ) (2.1) 

R 

Here Xr(U) is the group character of the element U in representation R, C 2 (R) is the eigen- 
value of the quadratic Casimir operator for the representation R, and A = g 2 (3L/4 where g 2 
is the gauge field coupling constant. The characters form an orthonormal basis on the vector 
space of irreducible representations of the gauge group and hence the propagator ( |2.1| ) is seen 
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u, 

Figure 1: The propagator between two open ends on the cylinder transports irreducible repre- 
sentations with quadratic Casimir dependent exponential damping. 



X.(U) X T (U) 




Figure 2: To insert matter into the theory the sewing property is frustrated by the introduction 
of a character of representation S at the sewing boundary. In the case of S being the trivial 
representation one recovers (2.3). 



to be a diagonal operator on this space. In the following we will consider the propagator as 
an operator, 

K[A] = e~ A ° 2 (2.2) 

The convolution of two propagators to form a single propagator follows from the sewing 
property (see fig.g), 

KlA^KlAz} = K[A X + A 2 ] (2.3) 

This relation follows directly from the orthogonality of group characters, which has the effect 
of multiplying the exponents of (|2.2|) in the naive way. This sewing is fundamental to using 
cylinders to build models of spherical or toroidal topology and to include external charges. 

The inclusion of a static charge of a particular irreducible representation S into the system 
|§ is accomplished by including character operator \s at a sewing site (see fig.0) . Here we will 
restrict ourselves to the case where all charges in the system are of the same representation 
but the extension for arbitrary configurations is straightforward. 

Introducing n, S representation loops each separated by a distance L\- ■ ■ L n , the operator 
of interest is 

n 

II [ ^ 2 ^ C2/4 Xs] (2.4) 
1=1 

If one identifies the ends of the cylinder and integrates over each Lj with an appropriate 
restriction, we arrive at an expression for the partition function of n charges on a circle with 
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Figure 3: Sewing n cylinders together with frustration and identifying the ends, leaves one 
with a torus corresponding to n charges on a circle. 



total circumference L = J2 L^. 

Z n = - [ L dLt . . . [ L dL N 8 (l - Tr f[ \ ^ L ^l\ s \ (2.5) 

n\Jo Jo V £i / t=x J 

It is worth noting here that this last expression also has the interpretation as a gas of Wilson 
loops on a space-time torus and with slight modification can be used as a starting point for 
investigating a gas of Wilson loops on the sphere. Taking the circumference of the circle, 
L to infinity effectively reduces the configuration space to an open line. Carrying out the 
now unrestricted integrations over {Li} we find that the partition function takes on a rather 
simple transfer matrix form 



1 n 
-7 Tr Il 



1 

Xs 



= ^-TrT n (2.6) 



Hence, the calculation of the thermodynamics of n static charges on the open line reduces to 
solving the eigenvalue problem for the operator T 

= = (2.7) 

For the purpose of finding thermodynamic quantities it is convenient to deal not with the 
constant volume ensemble as we have up to now, but rather a constant pressure ensemble. 
The change to a constant pressure ensemble can be carried out in a straightforward manner 
by introducing a pV term into the energy of the system with the result of shifting the energy 
per unit length. The resulting eigenvalue problem reads 

In the thermodynamical limit where n — > oo all information of the system is contained in the 
largest eigenvalue Ao of the operator T p . The remainder of this paper will involve finding Ao 
for the case of a gas of adjoint charges. 
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3 Low Density Limit: Group Theory 



The effective eigenvalue problem for the non-Abelian gas (|2.8|) was previously derived via a 
different approach by Nambu et al 0. As in that case (|2.8|) is equivalent to the following 
linear equation, 

= (aC 2 - qxs)^ = -P^ (3.1) 

which acts on the space of irreducible representations, so that ^ = J2r a RXR'i with a = g 2 /4 
and q = l/(/3A). The structure of this equation is the same as one would find in a quantum 
mechanics problem. The quadratic Casimir is a diagonal operator on the space of irreducible 
representations 

C 2 XR = C 2 (R) X R (3.2) 

which corresponds to the kinetic term. The role of the potential is played by the character xs 
which mixes the eigenvectors of the kinetic term. This can be easily seen by the multiplication 
rule 

Xr Xs = Xr®s = n rsXt (3.3) 

T 

Here N^ s is the fusion number which enumerates the occurrence of the irreducible represen- 
tation T in the Kronecker product of representations R and S. The only difference between 
quantum mechanics and the current situation is that we would like to solve for the eigenvalue 
of the transfer matrix problem A = l/(f3q) as a function of the pressure p as opposed to 
solving for the energy of the system as a function of the potential. 

As in the case of quantum mechanics one can begin to solve the eigenvalue problem by 
considering the symmetries of the system which will lead to conserved quantities. Here we are 
most interested in the symmetric properties of the 'Hamiltonian', H, under transformations 
which lie in the center of the gauge group. As we have seen, the presence of such a symmetry 
immediately leads to the phenomena of multiple vacua. The action of a transformation under 
the center of the gauge group is defined as 

Zxr = zrXr (3.4) 

Here zr is a representation of the center of the gauge group. Since for all compact Lie groups 
the center forms an Abelian subgroup, we can take zr to be a complex phase factor. The 
details of this phase factor depend on the structure of the center Z. For U(N), Z is isomorphic 
to U(l) hence Zr = e iaCl ^ where a G 1R and G\ is the integer valued linear Casimir operator. 
For the case of interest, SU(N), Z ~ TL^ and consequently zr = e 2mC ^( R )/ N _ This follows 
from the U(N) case with the restriction that Zr — 1. The question of whether H commutes 
with Z is reduced, since Casimir operators commute amongst themselves, to the question of 
whether 

[Xs,Z} = (l- e^ c ^ N ) X s = (3.5) 

The solutions to such a condition are clearly C\(S) = (mod N). In terms of irreducible rep- 
resentation S of the matter content of the theory this means that the full vacuum degeneracy 
is apparent only when C\(S) = (mod N). The simplest examples of such representations 
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are the trivial representation in which case the theory reduces to that of pure Yang-Mills in 
1 + 1 dimensions and the case where all matter is in the adjoint representation. These are the 
cases we will consider for the remainder of this paper. 

In the case of the adjoint gas where the center operator commutes with ( |3.1|) , in analogy 
with the conservation of eigenvalues of commuting operators in quantum mechanics we see 
that the eigenvalue of the linear Casimir is conserved (mod N) and is a good 'quantum number' 
(mod N). Consequently for the SU (N) adjoint gas there exists a family of N distinct solutions 
to the eigenvalue problem each of which we will label by k — 0, • - • , N — 1. For each value 
of k we have an isolated sector of the theory complete with a stable vacuum and an infinite 
tower of excited states. These are precisely the discrete '6*-vacua' of the model. 

In the limit where q — > 0, the eigenvalue problem in ( |3.1|) is reduced to that of the free 
1 + 1 dimensional Yang-Mills theory and we can easily identify the vacuum states of each 
sector. In this case the eigenvectors, \I/ of the transfer matrix, T p are simply the irreducible 
representations of the gauge group. Since we are interested in the case with adjoint SU(N) 
charges, the labeling of vacua introduced above with k = 0, • • • , N— 1 will be followed although 
strictly speaking the free theory has a countably infinite vacuum degeneracy. We denote 
by [k] the k th vacuum state which is the k(k — l)/2 dimensional completely antisymmetric 
fundamental representation of SU(N). Each of these fundamental representations is the 
lowest lying energy state of each of the k sectors and will serve as a starting point for a 
perturbative calculation of the eigenvalue problem for small q, or equivalently, small p in each 
sector. 

Having identified the ground states of each sector of the theory all we need to calculate, 
to lowest order in the pressure p, the solution of the eigenvalue problem (|3.1|) are the fusion 



numbers N]^ s ( |3.3|) . As we have seen, these are pure group theoretic quantities which detail 
the mixing effect of the potential on irreducible representations and, in particular, the anti- 
symmetric ground states [k]. With this information one can then proceed as in quantum 
mechanics and calculate with the unperturbed ground states the pressure p up to third order in 
q. Leaving the details of the evaluation of the fusion numbers and the perturbative expansion 
of the eigenvalue problem to the Appendix, we record the results for the SU (N) adjoint gas 
in the k th sector in Table 1. 

Now we would like to develop the equation of state for system. As is familiar from more 
physical gauge theories, the number of microscopic degrees of freedom n may not be the 
number of macroscopic degrees of freedom, n* . For example it is believed that in QCD 
pairs and triples of quarks are bound into observable mesons and bayrons, respectively. The 
equation of state per microscopic degree of freedom for the constant pressure ensemble is 

P<V> _ ^ ^ n* 

= p < v >= — T (3.6) 

n n 

where < V > is the expectation value of the total volume of the system which is canonically 
conjugate to p. It can be determined by inverting the relation p(q) and using the relationships 
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N > 2 
k = 


p — g + g + c^g J 

P = i-vf + o(p 2 ) 


N = 3 
k = 1 


^: ~ i 9 ~2 i 45 ~3 . /0/~4\ 
V = Q + 4<? + 16<? + °(<? ) 

p=l-fp+^ 2 + 0(p 3 ) 


N >3 
k = 1 


^ ^ i o N 2 z2 i /i JV 4 ^3 i r\i^\ 

v - <? + 2 N 2-iQ + 4 (7v2_i)2<? + <J{q ) 

p = 1 - 2^p + 4(^p 2 + 0(p 3 ) 


N > 3 
k = 2 


~ _ ^ i fiv i ^ i N ^n 2 i f n3 i ^ 2 i 2Af2 1 \F@ 1 

P — HT\2 N+l 1 7V-2^ 1 ^(/V-2)(jV+l) 1 (N-2) 2 1 (/V+l) 2 1 J« W \H ) 

1 /JV | JV | AT \ ~ 1 37V 4 +lCW 3 -13iV 2 -8iV-24 /\ r 2^2 , /Of™ 3 ^ 
P - 1 I 2 + JV+1 1 /V-2^ 1 4(JV-2) 2 (iV+l) 2 JV P + J 


N > 3 
k > 2 


| (N | TV | /V \~2 | / /V 2 | JV 2 | 2/V 2 2/V3 \~3 , Qfgt) 

n _l f/V | TV | TV N~ | / AT2 Af2 TV 2 5fc+7V+77Vfc+7V 2 N ~2 , 

P~ 1 v fc ^ 7V+1 1 N-k'r 1 H/V-fc) 2 (/V+l) 2 1 fc 2 (/V-fc)(/V+l) ^ ' 



Table 1: Tafr/e of the pressure p and the ratio of macroscopic to microscopic degrees of freedom 
p for the various vacuum states, k of the SU(N) adjoint gas. The primed variables are defined 
as x = -jr. Note the conjugation symmetry k —> N — k. 



between the thermodynamic variables 

^ dlogZ ^ dlogX d log q(p) 

< V >= -T — - = -nT— — = nT (3.7) 

op op op 

It is convenient to define p as the ratio of macroscopic to microscopic degrees of freedom 

n* d\ogq 

P = — = TTj 3.8 

n o log p 

Comparing with ( |3.6|) we see the fundamental importance of the quantity p. The dependence 
of p on p is tabulated for the different sectors of the SU(N) gas in Table 1. 

The results of these calculations deserve some comment. The most striking between the 
different sectors of the theory is the configuration of adjoint charges in the limit of vanishing 
pressure. For the k = sector we find that p = 1/2 — ■ • • and hence the adjoint charges in 
the system are bound pairwise in the low pressure limit. This behaviour is not surprising 
and is seen in both the U(l) Q and SU(N) || one dimensional (fundamental representation) 
Coulomb gases. What is different here in the k = sector is that the first corrections in 
pressure to this pair- wise binding come about with a negative sign and so the adjoint charges 
begin to form macroscopic configurations where the number of constituents is three of more. 
This is possible since adjoint charges are of course self-adjoint and an arbitrary number of 
them can form an observable charge singlet. 

When one moves to the cases when k > we see a distinct change in the vanishing pressure 
macroscopic structure of the theory. As explained previously for the Yang-Mills case 0, ||, 
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different sectors of the a 1 + 1 dimensional gauge theory are equivalent to considering a the 
theory with different constant background colour electric fields. For SU(N), each admissible 
background is given by one of the N fundamental representations which we label by the 
parameter k. In Table 1 we see that for a non-trivial background (k > 0) the adjoint charges 
of the system can interact with the background electric field and form stable, colour singlet 
configurations where they are the macroscopic degrees of freedom. In other words they act 
as free particles. 

4 High Density Limit: Fermions on the Circle 

The eigenvalue problem of (|3.1| ) can also be solved exactly in the limit of large values of q 
which corresponds to the limit of high pressure. This is most conveniently carried out by 
converting the group theoretic equation of (|3.1| ) to a linear differential equation with periodic 
coefficients. In this section the gauge group will be taken to be U(N) as this will simplify 
calculations. Recovering the results for SU(N) is a trivial step which will be noted at the 
appropriate point in the calculation. 

The starting point for converting the eigenvalue problem of (|3.1| ) to a differential equation 
is to consider the eigenvector $ as a linear combination of irreducible representations each 
labeled by N integers {rii, ■ ■ ■ ,tin}. These integers correspond to reduced row variables for 
the Young diagram associated with the irreducible representation, and satisfy the dominance 
condition 

oo > rii > n 2 > ■ ■ ■ > Un > — oo (4.1) 
The quadratic Casimir operator is diagonal in this basis with the action 

C2^,-,n N = \ (E n\ - N{N *~ 1} ) V„ lr ,» w (4-2) 

The action of the character xs on these states is, however, more complicated. This is partly 
due to the dominance restriction on the Young diagram. If one relaxes this restriction then 
the action of the character is simplified somewhat. In this case the action on a general state 
of a character in the adjoint representation is given by, 

N 

XAdflPni,---,n N = '^Pn 1 +S r ^ 1 -S s ^i,---,n N +8 r ^ N -8 a:N (4-3) 

r,s=l 

The result of this operation is to add unity to n r and then subtract unity from n s and sum 
over all r and s. Consequently the eigenvalue problem of ( |3 . 1| ) is a difficult recurrence type 
equation. This type of equation is most successfully dealt with by introducing the function 
ip{xi, ■ ■ ■ , xn), which is completely symmetric under the exchange of arguments, via a Fourier 
transform 

i lr , % = dxi--- dx N e l ^ n ^ijj(x 1 ,---,x N )A(x 1 ,---,x N ) (4.4) 
Jo Jo 
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The factor A({xj}) is a Vandermonde determinant defined as Yl^j^Xi — Xj) and as such is 
completely antisymmetric under the exchange of any two Xj's. This factor is included to force 
the integration to vanish identically for rii = rij, when i ^ j. In this way we can effectively 
impose the dominance condition ( |4.1| ) of the reduced Young diagram variables in the Fourier 
domain. 



Acting with the Casimir and the adjoint character on ( |4.4| ) we find the transfer matrix 
problem of ( |3.1| ) is equivalent to the second order linear partial differential equation with 
periodic coefficients 



1 N d 2 N 
--T — - — (N 2 



1)-* 

a 



N + 2 2J cos {xi — Xj 



i<j 



Aip = --Aip (4.5) 

a 



In this form some of the features of the adjoint non-Abelian Coulomb gas are more apparent. 
For example, the center operator Z has a simple interpretation in the Fourier domain 



Z 



iaCi 





N 




\ 4r d 1 


exp 


ia^n k 
. fe=i 


= exp 


a ^-dx h 

. fc=l ux k_ 



(4.6) 



This is exactly the same structure as the translation operator in quantum mechanics. Here 
Z generates uniform shifts of the coordinates {x{\ — > {xi + a}. It is easy to verify that (|4.5|) 
has this symmetry and from the arguments of the previous section we expect a continuum of 
vacua for the adjoint gas with U (N) gauge group. Additionally with periodic coefficients and 
the restriction to completely anti-symmetrized wavef unctions, Aip, the eigenvalue problem is 
equivalent to that of non-relativistic fermions in a periodic potential. This correspondence is 
familiar from matrix models [H| and is exploited in the solution of the large non-Abelian 
Coulomb gas [JTTJ, [12], [13 |. 

In the limit of high densities, or equivalently high pressure, the wave functions are localized 
about the minima of the potential. Expanding the potential about the local minimum at 
{xi — xj = 0} leads to the coupled harmonic oscillator problem 



N 



d 



N 



oEi-^-i) + !E(^- 



a 



N 2 q 



a 



Aip 



Performing the change of variables to the orthonormal basis {ui} given by 
((JV -n) x n - {x n+ i + x n+2 + ■■■ + x N )) 



—Aip 
a 



(4.7) 



u n 



n 



1...JV-1 



1 



'(N -n)(N -n + 1) 
(xi + x 2 + ■ ■ ■ + x N ) 



(4.8) 



diagonalizes the problem. In this basis the decoupled oscillator problem is 

1 



N Q2 



2 § du 2 ' 2 V a 



' 8=1 



ip = Eip 



(4.9) 
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where, 

E = " 1) - P ~ N2q (4-10) 

At this point we note that there is no potential for the coordinate which is to be expected 
as it describes a center of mass coordinate in the change of variables ( |4.8|) and the original 
potential in ( |4.5| ) depends only on relative not absolute positions. Consequently the un 
dependence of the problem is only through a phase. In the case of SU(N) this center of 
mass coordinate is restricted but otherwise behaves exactly as in U(N), entering only as 
a phase. Regardless of the details of this phase we will see it does not contribute to the 
high pressure equation of state of the adjoint gas. The other modes corresponding to the 
coordinates {ui ■ ■ -un-i} have degenerate frequencies which are easily read off the diagonal 
form 

u N = J— (4.11) 
V a 

Knowing the normal modes of the eigenvalue problem we are now in a position to find the 
ground state solution of (|4.5|) which will correspond to the dominant eigenvalue of the transfer 
matrix problem. A solution which satisfies the requirement of antisymmetry with respect to 
permutation of the coordinates is given by, 

Ai> ~ Y[(xi - Xj) exp ji-^5>ijexp j-i^^fo -x,) 2 j (4.12) 

The parameter M is an integer associated with the center of mass coordinate and contributes 
a constant to the energy eigenvalue. Notice that the potential exponentiated, this is a direct 
consequence of the normal modes having degenerate energy. The above state has iV 2 — 1 
quanta of energy, 

E = (Ar2_!)^ + I M 2 ( 413 ) 
so that the pressure is given, up to an irrelevant constant by, 



l = ^.^-l)M (4,4) 
a a V 2a 

Inverting this relation to find q(p) and using the definition fl3.8|) we find, to leading order, the 
ratio of macroscopic to microscopic degrees of freedom 

Consequently, in this high pressure limit the adjoint charges are the macroscopic degrees of 
freedom and, as in the k > low pressure cases, can be interpreted as being free. The most 
striking feature of (|4.15|) is the absence of any k dependence. This follows from the fact 
that the information of the center of mass coordinate appears only as a phase contributing 
the additive constant M 2 to the eigenvalue problem which is negligible in the limit of large 
pressure. In terms of physics the high pressure adjoint gas effectively screens all colour electric 
fields over large distances and so the fundamental colour electric fields associated with the 
different vacuum sectors are washed out by the adjoint degrees of freedom. 
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5 Discussion and Conclusions 



For any finite N the SU(N) non-Abelian Coulomb gas does not admit phase transitions, 



however, in the limit of infinite N it has been shown [11], [L2], |T3| that a phase transition 
develops for the trivial vacuum sector. In those works, it was found that the adjoint 'quarks' 
remain bound strictly pairwise, to form 'mesons', up to some critical pressure at which a 
first order phase transition occurs and the mesons disassociate into free quarks. If one takes 
the infinite N limit [] in the trivial vacuum sector (k = 0), for small pressures the ratio of 
macroscopic to microscopic degrees of freedom is (see table [I|), 

P=\-^l + 0{P) (5.1) 

which appears to contradict the prediction of the large iV calculations. Equation fl5.1| ) predicts 
that the quarks are bound pairwise only at zero pressure, and can form bound structures 
with more components slightly away from zero pressure. This discrepancy can be explained 
by noting that in the infinite N limit all pressures are naturally measured in units of iV 2 . 
Consequently the statement that quarks are bound pairwise from zero pressure up to some 
critical pressure is misleading. The correct statement is: for pressures of order N 2 but less 
than the critical pressure, the quarks are bound strictly pairwise. In this large N limit nothing 
can be said about pressures of order unity, where our current computations are valid. We see 
then that the infinite N computations of [IT , 12, [H| are complementary to our own analysis. 



To state it another way, our computations cannot say anything about the affect that multiple 
(vacuum) sectors in the theory have on the phase transition at large N. However in the limit 
of infinite N our calculations are perfectly valid in the low pressure region. Taking the infinite 
N limit of the density for N > 3, k > 2 one finds (see table P, 



where k = ON . This demonstrates the explicit dependence on vacuum sector in the small 
pressure region. 

In the high pressure limit the picture is much simpler. The infinite N limit of the exact 
expression for the ratio of macroscopic to microscopic degrees of freedom p is given by 

l + (l + 8p)-5 

P = 4p — - (5.3) 

1 + 4p+ (l + 8p)3 

Here we have taken the pressure to be of order N 2 by setting p = p N 2 before taking the 
infinite N limit. This ensures that the expression for the density is in the correct pressure 
scale. Taking the pressure to be of order N 2 allows us to probe near the phase transition 
region, but not the transition itself, since our high pressure approximation of the last section 
has neglected the topology of the problem which drives this transition. For instance we find 



3 In this limit the product aN is kept at order unity, and we use p to denote the reduced pressure 



N 
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that at a pressure of p — | the density of the gas is exactly one half. This gives a lower 
bound on p cr it. since the phase transition must occur before the gas confines. Also at very 
large pressures we have the prediction 



P=l-^= + -\ (5.4) 




which agrees with the large N computations. 

Finally, some comments on the application of our results to the case of two dimensional 
U(N) Yang-Mills theory coupled to adjoint matter. For adjoint representation fields with 



mass m the high temperature effective potential can be calculated |14|, [15], [L6] 



V(xi — Xj) ~ — - — K\(m{5k) cos [k(xi — Xj)} (5.5) 
fc=i k 

Here K\ is a modified Bessel function. In the limit of large mass such that mj3 is large it can 
be shown that the effective potential ( |5.5| ) reduces to the form of the transfer matrix problem 
(|4.5| ) we have considered previously. For the case of finite mass, the full high temperature 
potential is difficult to deal with but it clearly has a symmetry under translations of the 
form {xi} — > {xi + a}. As we have seen previously, such a symmetry is synonymous with 
invariance of the model under center transformations ( |4.6|) . This strongly suggests that our 
results for the classical adjoint gas carry forward to massive adjoint two dimensional QCD at 
high temperature and we should expect multiple vacua to exist in that theory. Moreover, the 
physics in each sector of the theory will depend on the superselection parameter which labels 
the sectors. 

In conclusion, from consideration of the thermodynamics of a system of static adjoint 
representation charges interacting via SU(N) colour electric fields in 1 + 1 dimension we have 
shown that the physics depends on the discrete vacuum index k. We have solved the model 
in two regions: low and high pressure. In the low pressure regime, which is equivalent to low 
particle density, the constant pressure equation of state was shown to have strong dependence 
on k. In the limit of high pressures/particle densities the dependence on k was shown to 
become trivial, and does not enter into the equation of state for the adjoint gas. This is 
attributed to the screening nature of the high pressure limit, which washes out any global 
structure like a vacuum index. 
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Appendix 



Here we present the details of the perturbative calculation of the lowest lying eigenvalue of 
equation (|3.1| )for the adjoint gas. This can be accomplished via the familiar time independent 
perturbation theory of quantum mechanics. The problem we are faced with is 

(aC 2 - qXAd)^ = -p9 (A.l) 

As noted previously the eigenvectors of the Casimir operator C 2 are the irreducible repre- 
sentations of the gauge group so a convenient basis in the limit as q — > is \I/ = YJ clrXr- 
Also we have seen that the ground state of each sector of the theory is given by one of the N 
fundamental representations, {[k]} of SU(N). Consequently our objective is to calculate the 
perturbations to the pressure p for small q for ( |A. 1| ) about each ground state [k]. In order to 
carry out this calculation we need to know the matrix elements of the potential in the basis 
of the irreducible representations. This information follows from the Kronecker product of 
the adjoint representation with our chosen basis 

XAdXR = XAd®R = J2 N Ad RXT (A.2) 
T 

Here Yj d R is the fusion number enumerating the occurrence of the irreducible representation 
T in the product of the adjoint representation Ad and R. As in quantum mechanics we can 
easily calculate the corrections to p up to third order in q using the unperturbed basis of 
irreducible representations. 



PW ~ ^ Ad[k]+ q ^ cm _ C2{R) 



+q 3 



ly Ad S ly Ad R ly Ad [k] /» o\ 

b&W (C 2 ([k\) - C 2 (R)) (C 2 ([k]) - C 2 (S)) [ ■ } 



, N [k] sr Ad w 

Ad[ \% ] \C 2 ([k])-C 2 (R)_ 



+ 0(g 4 ) 



It should be noted that we have left out a constant, sector dependent, background contribution 
to the pressure. 

In order to give the details of the calculation we need a notation to label the irreducible 
representations. The one we will use here is given by the column variables [7711,772,2,...] 
of the Young diagram associated with the representation. For example the antisymmetric 
combination of k, N dimensional fundamental representations in SU(N) - the ground state 
of the k th sector - corresponds to a Young diagram with a single column of k boxes: [k]. 
Another example which appears in all calculations is that of the adjoint representation (Ad) 
which in column variables is given by: [N — 1, 1]. In this notation the quadratic Casimir for 
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a representation [m 1; m 2 , . . .] is given by 



C 2 (R) = V 



( 1-F 
1-1 
1-f 



1-f 
2(1- f) 
2(1- f) 



1-f 
2(1 -f] 
3(1-10 



Vi 



TV-1 



2(1 



JV-l - 



3(1 



TV-1 ' 
TV - 



1 - 

2(1- 
3(1- 



TV-1 
TV 

TV-1 ' 
TV - 
TV-1 ' 

TV - 



(tV-1)(1 



TV-l 
TV 



V T --(N 2 
12 v 



)/ 



(A.4) 



where 



V = (mi + 1, m 2 + 1, . . . , mjv-i + 1) 



The remaining task is to compute the relevant fusion numbers. We begin by presenting the 
results of the calculations in Table 2. Here we record the fusion numbers N^ d R = N^ d s for 
the representations R and 5* of importance in the calculation (|A.3|) of the pressure. Each sub- 
table corresponds to a different background k for SU (N) since the details of the calculation 
of fusion numbers in general depends on k and N. These results are only good for k < N/2 
where the remainder of the cases can be found via the symmetry of the eigenvalue problem 
under conjugation k — ► tV — k. For completeness we present the details of the second table 
for k = 1 and tV = 3. This is the familiar case of SU(3) and via common tensor or Young 
diagram methods the Kronecker products of the 8-dimensional adjoint representation (Ad) 
with the lowest lying representations can be calculated. In dimension notation we have 



®3 = 3©6©15 
®6 = 3©6©15©-- 
g) 15 = 3©6©2x 15 



(A.5) 



In the last two products we have ignored higher representations which do not contribute to 
0(q 3 ) in ( |A.3| ). Converting to our column notation 

(A.6) 
(A.7) 



3 = 


[1] 


6 = 


[2,2] 


8 = 


[2,1] 


15 = 


[2, 1, 1] 



we have the results of the second sub-table in Table 2. 
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k — 0, TV > 2 


[0] [TV -1,1] 




[0] 

\M 111 


1 

1 2 






k — 1, TV = 3 


[1] [2,2] [2,1,1] 




[1] 
[2,2] 

[2, 1, 1] 


1 1 1 
1 1 1 
1 1 2 






k — 1, TV > 4 


[1] [TV -1,2] [TV -1,1,1] 




[1] 

[TV -1,2] 
[TV -1,1,1] 


1 1 1 
1 2 1 
1 1 2 






fc = 2, TV > 3 


[2] [TV -1,3] [1,1] [TV -1,2,1] 




[2] 

[TV -1,3] 

[1,1] 
[TV -1,2,1] 


111 1 
12 1 
10 1 1 
111 3 


k > 3, TV > 3 


[k] [N-l,k + l] [A; -1,1] [TV-l,/c,l] 


[k] 

[TV - 1, fc + 1] 
[A; -1,1] 
[TV -l,k,l] 


11 11 

12 1 

10 2 1 

11 13 



Table 2: Table of relevant fusion numbers N^ d R = N^ d s for the calculation of the pressure 
of the adjoint gas. Note these results hold only for k < TV/2 with the other cases given by the 
symmetry k — > TV — k. 
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